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ABSTRACT 

We measure the faint end slope of the galaxy luminosity function (LF) for cluster galaxies at 
1 < z < 1.5 using Spitzer IRAC data. We investigate whether this slope, a, differs from that of 
the field LF at these redshifts, and with the cluster LF at low redshifts. The latter is of particular 
interest as low-luminosity galaxies are expected to undergo significant evolution. We use seven high- 
rcdshift spectroscopically confirmed galaxy clusters drawn from the IRAC Shallow Cluster Survey 
to measure the cluster galaxy LF down to depths of M* + 3 (3.6/im) and M* + 2.5 (4.5/im). The 
summed LF at our median cluster redshift (z = 1.35) is well fit by a Schechter (1976) distribution 
with a^Sfj-m — —0.97 ± 0.14 and a4.5 Mm = —0.91 ± 0.28, consistent with a flat faint end slope and is 
in agreement with measurements of the field LF in similar bands at these redshifts. A comparison to 
a in low-redshift clusters finds no statistically significant evidence of evolution. Combined with past 
studies which show that M* is passively evolving out to z ~ 1.3, this means that the shape of the 
cluster LF is largely in place by z ~ 1.3. This suggests that the processes that govern the build up 
of the mass of low-mass cluster galaxies have no net effect on the faint end slope of the cluster LF at 
z<1.3. 

Subject headings: galaxies: clusters: general, galaxies: evolution, galaxies: formation, galaxies: lumi- 
nosity function 



1. INTRODUCTION 

Many studies have shown that the low-mass cluster- 
galaxy population evolves substantially at low redshift. 
For instance, Cowie et al. (1996) first recognized that 
star formation happens primarily in high-mass systems 
at high redshift and low-mass systems at low redshift, 
a fact which has been studied extensively since (see, for 
example, Panter et al. 2007; Mobasher et al. 2009; Chen 
et al. 2009; Villar et al. 2011). Moreover, it is well known 
that cluster galaxies undergo morphological transforma- 
tion at low redshift, with many cluster members trans- 
forming to lenticular galaxies at low redshift (Dressier 
et al. 1997; Desai et al. 2007; Wilman et al. 2009). There 
is also substantial evidence that the low-luminosity red- 
sequence galaxy population grows substantially in clus- 
ters since at least z ~ 1 (Stott et al. 2007; Lu et al. 2009; 
Rudnick et al. 2009; Lemaux et al. 2012). 

Taken together, these facts demonstrate that the low- 
mass cluster galaxies are actively evolving and forming 
since z ~ 1. Therefore, by comparing low mass, z = 
galaxies with their high-redshift progenitors we can po- 
tentially constrain the processes important in galaxy for- 
mation and evolution. This can be done by studying 
individual galaxies (through their star formation rates, 
stellar masses, morphological types, and structural prop- 
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erties) or by studying galaxy populations (through their 
luminosity and mass functions). 

In particular, the near-infrared luminosity function 
(NIR LF) can be used to study the stellar mass growth 
of a galaxy population, as the rest-frame NIR is a good 
proxy for stellar mass (Muzzin et al. 2008). In clusters, 
the NIR LF has been used extensively to study the as- 
sembly of the most massive cluster galaxies. Such studies 
have found that the massive end of the NIR LF evolves 
passively out to z ~ 1.3, suggesting that the bulk of the 
stellar mass of these galaxies is in place at high redshift 
(Andreon 2006; Strazzullo et al. 2006; De Propris et al. 
2007; Muzzin et al. 2008; Mancone et al. 2010). In ad- 
dition, in Mancone et al. (2010) we found statistically 
significant deviations from passive evolution at z > 1.3 
which we could only explain with ongoing stellar mass 
assembly at these redshifts. 

Most attempts to probe the faint end of the cluster lu- 
minosity function (LF) at high redshift have been limited 
to studying the red sequence. Such studies have found a 
deficit of faint and red cluster members at high redshift 
when compared to their low-redshift counterparts (De 
Lucia et al. 2004; Stott et al. 2007; Lu et al. 2009; Rud- 
nick et al. 2009; Lemaux et al. 2012). This could mean 
that low-mass cluster galaxies undergo substantial mass 
growth at low redshift, or simply that low-mass cluster 
galaxies are still blue at high redshift and have not fin- 
ished transitioning onto the red sequence (Lemaux et al. 
2012). Differentiating between these two cases requires 
measuring the LF of all faint cluster members. Previ- 
ously, Strazzullo et al. (2010) was the only study to do 
this, finding a faint-end slope consistent with fiat. How- 
ever, they did not compare their results to low-redshift 
clusters to determine the implications for the stellar mass 
growth of low-mass cluster galaxies. 

In this paper we measure the 3.6 and 4.5 /jm LF of high 
redshift (1 < z < 1.5) galaxy clusters. Our measure- 
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TABLE 1 

Cluster Member Summary 



Cluster RA Dec z # Members 

ISCS J1432.4+3332 14:32:29.18 33:32:36.0 1.112 26 

ISCS J1434.5+3427 14:34:30.44 34:27:12.3 1.238 19 

ISCS J1429.3+3437 14:29:18.51 34:37:25.8 1.261 18 

ISCS J1432.6+3436 14:32:38.38 34:36:49.0 1.351 12 

ISCS J1433.8+3325 14:33:51.13 33:25:51.1 1.369 6 

ISCS J1434.7+3519 14:34:46.33 35:19:33.5 1.374 10 

ISCS J1438. 1+3414 14:38:08.71 34:14:19.2 1.414 16 



ments trace the rest-frame NIR, where the LF is known 
to correlate well with stellar mass. Most importantly, 
our data are deep enough to constrain a, the faint-end 
slope of the LF. This, combined with low-redshift results 
from the literature, allows us to measure the stellar mass 
buildup of the low-mass cluster galaxy population over 
a redshift range when these galaxies are known to be 
actively evolving. 

This paper is structured as follows. Section 2 describes 
our data. Section 3 presents our method for measuring 
the galaxy cluster luminosity function and gives our re- 
sults. In Section 4 we compare our results to low-redshift 
clusters and the field. Our conclusions are presented in 
Section 5. All magnitudes are on the Vega system, and 
we assume a WMAP 7 cosmology (Komatsu et al. 2011; 
O m = 0.272,17a = 0.728, h = 0.704) throughout. All 
SPS model predictions are generated using EzGal (Man- 
cone & Gonzalez 2012). 

2. DATA 
2.1. Cluster Sample 

The clusters from this study are part of the IRAC Shal- 
low Cluster Survey (ISCS) (Stanford et al. 2005; Elston 
et al. 2006; Brodwin et al. 2006; Eisenhardt et al. 2008), 
a catalog of clusters identified as 3-D overdensities using 
photometric redshifts in the 8.5 deg 2 Bootes field. Fur- 
ther work with the high-rcdshift (1 < z < 1.5) clusters 
in the ISCS has included deep (1000s) IRAC imaging, 
spectroscopic followup, and HST imaging. Seven of the 
spectroscopically confirmed, high-redshift ISCS clusters 
have both deep IRAC imaging and ACS F775W imaging. 
It is this subsample of ISCS clusters that we use to study 
the LF of high-redshift galaxy clusters. We supplement 
the 1000 seconds of targeted IRAC observations for each 
cluster with imaging from the Spitzer Deep, Wide-Field 
Survey (SDWFS, Ashby et al. 2009) which has a me- 
dian exposure time of 420 seconds throughout the Bootes 
field. This gives a total observing time of roughly 1400s 
per cluster in all four IRAC bands. 

We list in Table 1 our clusters along with their posi- 
tions, redshifts and number of spectroscopic members. 
ISCS J1438. 1+3414 has a published X-ray mass esti- 
mate of log(M 2 L o ^/M ) = 14.35^014 which comes from a 
143 ks Chandra exposure (Andrcon et al. 2011; Brodwin 
et al. 2011). All of these clusters (with the exception 
of ISCS J1433. 8+3325) have a weak-lensing mass esti- 
mate from Jee et al. (2011), with masses in the range of 
14.40 < log(M^ L /M Q ) < 14.73. 

2.2. Comparison Fields 



Given the limited spectroscopic redshifts in these fields, 
statistical background subtraction is required to recover 
the underlying LF. A statistical background subtraction 
involves measuring the number counts of galaxies in a 
field region and subtracting it from the number counts 
of galaxies near the cluster. This technique has been used 
successfully by de Propris et al. (1998), Lin et al. (2004), 
Muzzin et al. (2008), and Mancone et al. (2010). This 
requires a survey with IRAC imaging of at least the same 
depth as our cluster images as well as ACS F775W imag- 
ing. For this purpose we select the GOODS North and 
South (Dickinson et al. 2003) fields. We downloaded the 
latest fully reduced 3.6 and 4.5/im Spitzer IRAC images 
taken of the GOODS fields. We also retrieved the lat- 
est HST ACS F775W catalogs from the GOODS survey 
Throughout this paper we refer to the GOODS fields as 
our control fields. 

2.3. Data Reduction and Processing 

We produced IRAC mosaics of all seven clusters 
by combining data from our own programs (PID78, 
PID30950) with that from SDFWS, following procedures 
identical to those described in Ashby et al. (2009). This 
included the manner in which outliers were rejected and 
in the way the individual IRAC frames were prepared for 
mosaicing by first removing the residual images arising 
from earlier exposure to bright sources. We generated 
catalogs by running our fully reduced 3.6 and 4.5^m 
IRAC images (for the clusters and the control fields) 
through Source Extractor (Bcrtin & Arnouts 1996) in 
single-image mode. We used All diameter aperture mags 
that were aperture-corrected to total mags by comparing 
All and 24// diameter aperture magnitudes for bright, un- 
saturated stars in our images. We used stars to measure 
the aperture corrections because galaxies at these red- 
shifts are typically unresolved in IRAC imaging. This 
gave aperture corrections of —0.32 (—0.34) magnitudes 
in 3.6 (4.5)/^m for our cluster images and —0.31 (—0.32) 
mags for our control images. For reference, the difference 
between the aperture corrections of the cluster and con- 
trol fields is smaller than the uncertainty of the absolute 
flux calibration for IRAC images (Reach et al. 2005). To 
verify our calculated aperture corrections we compared 
our All aperture corrected magnitudes to the All aper- 
ture corrected magnitudes from SDWFS, and found very 
small systematic offsets (< 0.03 mags). 

The ACS imaging for our clusters was obtained as part 
of the HST Cluster Supernova Survey, and the reduction 
of the images is described in detail in Suzuki et al. (2012). 
We ran the reduced ACS F775W images through Source 
Extractor and used MAG_AUTO to calculate the F775W 
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Fig. 1. — Measured 50% completeness limits for our images in 3.6 (left) and 4.5^m (right). In both panels the solid histograms show the 
completeness limits for our cluster fields, and open histograms show completeness for our control fields. The vertical dashed line denotes 
the 50% completeness limits from the SDWFS survey (Ashby et al. 2009). 



magnitude of our galaxies. We used Our F775W pho- 
tometry to perform an optical— NIR color cut to remove 
contaminants from our LF (Section 3.1). For our control 
we used the ACS F775W MAG.AUTO values from the 
GOODS catalogs, which were also generated with Source 
Extractor. 

Next we calculated completeness as a function of mag- 
nitude at 3.6 and 4.5^m for each cluster image and 
each control image separately. We approximated our 
galaxies as point sources due to the coarse IRAC point 
spread function (PSF). We generated 24,000 artificial 
point sources for each image, uniformly distributed be- 
tween 13 and 25 mags. Our artificial point sources were 
simply copies of the PSF for each image, which we gen- 
erated by median combining unsaturated stars taken di- 
rectly from each image. We added these sources to the 
original images ten at a time, ran Source Extractor again 
for each new image, and finally calculated the recovery 
rate as a function of magnitude for a given image and 
filter. 

Figure 1 shows a histogram of the measured 50% point 
source completeness limits for each of our cluster images 
and control images at 3.6 (left) and 4.5^im (right). For 
comparison the vertical black line denotes the 50% com- 
pleteness limit in each band from the SDWFS survey. We 
only want to fit for the cluster galaxy LF when the com- 
pleteness of all galaxies in all clusters is at least 50%. 
Therefore we limit our LF fitting procedure to the the 
brightest 50% point source completeness limit for all our 
clusters, which is 20.37 (19.60) mags in 3.6 (4.5)^m. The 
50% completeness limit for our cluster images is ~0.75 
mags fainter than for the SDWFS images and ~1 mag 
brighter than our control images. 

3. OBSERVED LUMINOSITY FUNCTION 

3.1. Optical-NIR Color Cut 

We use a simple color cut to remove stars and low- 
rcdshift galaxies from our catalogs and increase the 
signal-to-noise ratio of our high-redshift cluster galax- 
ies. Our color cut is designed to include the blue cloud, 
as excluding part of it would induce a systematic bias in 



our measurement of the faint-end slope. We choose our 
color cut by using the Bruzual & Chariot (2003) stellar 
models to create a model of a star-forming galaxy with a 
color on the blue side of the blue cloud, consistent with 
Lemaux et al. (2012). We then use this same model to es- 
timate the color of the bluest star forming galaxies in our 
clusters, finding F775W-[3.6] > 3.5 and F775W-[4.5] > 
3.75. We use these values for our color cut, and note that 
our final results are not sensitive to our exact choice be- 
cause our results change by less than our random errors 
for a wide range of color cuts (2.5 < F775W— [3.6] < 5). 

Stellar contamination is a potential issue for our sam- 
ple as our cluster and control fields are at different galac- 
tic latitudes. A [3.6]- [4.5] cut could effectively remove 
stars from our sample, but is not possible because there 
is little overlap between the 3.6 and 4.5 /an imaging of 
our control regions. Limiting our sample to areas with 
3.6 and 4.5 /xm imaging would remove about 75% of our 
control region. However, stellar contamination is effec- 
tively removed via our optical-NIR color cut. To verify 
that stellar contamination is not an issue for our results 
we calculate the expected colors for local stellar popula- 
tions. To do this we download stellar isochrones using 
the CMD 2.3 software 8 which includes the latest stel- 
lar modeling details from a number of sources (Bonatto 
et al. 2004; Girardi et al. 2008; Marigo et al. 2008; Gi- 
rardi et al. 2010). For a low metallicity (Z = 0.008) 
model, which is relevant to the Galactic halo, no star of 
any age has F775W— [3.6] > 4. Solar metallicity stars 
have F775W— [3.6] ~ 5, at the reddest. Therefore, only 
the tip of the RGB and AGB for old stars extend red- 
der than the color cut. As such, only a small fraction 
of stars might remain after the cut, meaning that stel- 
lar contamination is not an issue. This is confirmed by 
the fact that our results do not change even when using 
color cuts as red as F775W— [4.5] = 5, which removes 
all stellar contamination. We additionally remove from 
our catalogs all objects with CLASS.STAR > 0.8 in the 
ACS catalogs. We find that this has a negligible impact 
on our results, again showing that stellar contamination 

8 http://stev.oapd.inaf.it/cmd 
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[3.6] [4.5] 

Fig. 2. — Binned and background-subtracted luminosity functions for 3.6 (left) and 4.5^tm (right). Median cluster redshift and best fitting 
Schcchtcr parameters are displayed in the top left. The solid curve shows the best fit. The dashed vertical line illustrates the magnitude 
limit for the fit, which is determined by the 50% completeness limit of the shallowest cluster image. 



is not an issue for our sample. 

3.2. LF Fitting Procedure 

We design our methodology so that we can perform an 
unbinned fit to the cluster member LF, and so that the 
background subtraction is done in a way equivalent to a 
subtraction in observed space. We generate an individual 
cluster LF and control LF for every cluster. The clus- 
ter LF is simply composed of the galaxies in the cluster 
image which passed our various cuts (Section 3.1), are 
within 1.5 Mpc of the cluster center, and are outside 
of the heavily blended cluster cores (typically ~100kpc). 
The latter restriction also removes the BCGs from our 
sample, which are known to not follow a Schechter dis- 
tribution. For each cluster we use a Bruzual & Char- 
lot (2003) SPS model to calculate the fc-correction and 
distance modulus correction needed to move a passively 
evolving galaxy (zf — 3, Chabrier IMF, solar metallic- 
ity) from the cluster redshift to the median redshift of our 
cluster sample (£=1.35). We then apply this fc-correction 
and distance modulus correction to all galaxies in the 
cluster LF. Next we build a control LF for each cluster 
in a similar fashion. We select all galaxies in the con- 
trol images which pass our cuts, weight them according 
to the relative area of the cluster and field images, and 
apply the exact same fc-correction and distance modulus 
correction that we applied to the cluster LF to all the 
galaxies in the control LF. We do this so that the same 
transformation has been applied in the same way to the 
cluster and control galaxies, and therefore when we sub- 
tract the control LF from the cluster LF the subtraction 
is effectively done in observed space. 

This procedure gives us an unbinned cluster and con- 
trol LF for each cluster. We then combine the individ- 
ual cluster and control LFs into a composite cluster and 
composite control LF, which we use to measure the LF of 
cluster members. We parameterize the luminosity func- 
tion of cluster members as a Schechter (1976) luminosity 
function and measure the best fitting Schechter parame- 
ters with maximum likelihood fitting, similar to the pro- 
cedure used in Mancone et al. (2010). This procedure 
requires an analytical representation for the contribution 
from the control region so we bin our composite control 
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Fig. 3. — Confidence regions for the Schechter fits to our 3.6/im 
(filled) and 4.5/iin (dashed) cluster LFs. Contours represent the 
1, 2, and 3a confidence regions in a vs M* space. Filled circles 
denote the best fit Schechter parameters 

LF by magnitude, correct for photometric incomplete- 
ness, and fit a third order polynomial to it in log space. 
We then use maximum likelihood fitting to fit the sum 
of a Schechter luminosity function (the cluster member 
LF) and the fitted composite control LF to the composite 
cluster LF. We use a downhill simplex algorithm (Press 
et al. 2007) to maximize the likelihood as a function of 
$*, to*, and a, and fit all galaxies brighter than the 50% 
completeness limit for the clusters (Section 2.3). 

3.3. Results 

Figure 2 shows the control-subtracted cluster LF and 
the Schechter fit to the cluster member LF for 3.6/im 
(left) and 4.5/Km (right). Maximum likelihood fitting 
gives a fit to the LF without binning, but for plotting pur- 
poses we show the binned and control-subtracted cluster 
LF in Figure 2, which is the binned difference between 
the composite cluster LF and the composite control LF. 
In each panel the solid curve shows the best fit while the 
dashed vertical line illustrates the magnitude limit used 
for the fit. Figure 3 shows the 1, 2, and 3<7 contours in 
M* vs. a space derived from our measured likelihoods. 
We also report the binned LF values in Table 3, although 
we note that our fit was to the unbinned data. 
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TABLE 2 

Best Fitting Schechter Parameters 



(z) # Clusters M* ^ m 


C*3.6nm 


M 4*.5 Mm 


«4.5(im 


1.35 7 17.44 ±0.30 


-0.97 ±0.14 


17.10 ±0.43 


-0.91 ± 0.28 



TABLE 3 
Binned LFs 



Mag 



# [3.6] 



# [4.5] 



15.25 
15.75 
16.25 
16.75 
17.25 
17.75 
18.25 
18.75 
19.25 
19.75 
20.25 



-1.32 ±0.50 
10.64 ±9.10 
0.87 ± 13.97 
32.86 ± 23.12 
95.61 ±32.39 
128.29 ± 37.82 
136.13 ±41.68 
181.98 ± 48.57 
145.11 ± 51.32 
270.40 ± 62.50 
253.25 ±71.37 



7.09 ± 7.02 
-3.31 ± 10.12 
24.34 ± 20.76 
88.94 ±30.76 
114.71 ±37.00 
118.11 ±41.72 
143.18 ±45.37 
214.39 ± 54.09 
182.13 ±59.06 



We also measure uncertainties using bootstrap resam- 
pling, as such an error estimate is more sensitive to sys- 
tematic uncertainties caused by cluster-to-cluster varia- 
tions. We generate realizations of the LF by randomly 
selecting seven clusters from our sample and repeating 
our LF fitting procedure with the new cluster sample. 
The cluster selection is done with replacement, which 
means that an individual cluster can be selected more 
than once when generating a new cluster sample. This 
allows us to probe any systematic uncertainty caused by 
cluster-to-clustcr variation, as this process effectively ap- 
plies random weights to the clusters while fitting. We 
perform 100 realizations of the cluster LF and take the 
standard deviation of the fitted M* and a parameters 
as our random uncertainties. Our best fitting Schechter 
parameters and bootstrap errors are listed in Table 2. 
We note that our measured bootstrap uncertainties agree 
well with the contours in Figure 3. 

The relatively small sizes of the GOODS fields means 
that cosmic variance in our control fields could be an 
additional source of systematic uncertainty. To verify 
that cosmic variance is not strongly biasing our results 
we redo our fit but use just one of the GOODS fields 
for our control sample, and then redo it again using the 
other. In each case the best fitting value of a changes by 
~0.1 in both filters, which is smaller than our measured 
errors. Therefore, cosmic variance is unlikely to be a 
dominant source of uncertainty in our results. 

4. DISCUSSION 

4.1. High-Redshift Comparison 

For a basic consistency check we compare to our re- 
sults from Mancone et al. (2010). In Mancone et al. 
(2010) we measured M| 6Atm and M| 5 out to z = 1.8 
using the slightly shallower SDWFS data and a statis- 
tical background subtraction. As such the methodology 
is very similar, the filters are the same, and the seven 
clusters studied herein were also included in Mancone 
et al. (2010). Due to our shallower data in Mancone 



et al. (2010) we fixed a and reported fitted M* values for 
a = —0.6, —0.8, and —1.0 in redshift bins from z = 0.3 to 
z = 2.0. The median redshift of the clusters in this study 
(z = 1.35) fall directly between the z = 1.24 and z = 1.46 
bins from Mancone et al. (2010). Therefore we compare 
our fitted M* values to the average of the M* values for 
the z = 1.24 and 1.46 bins with a = —1.0, which gives 
M 3 * 6jum = 17.42 ±0.1 and M 4 * 5/im = 16.92 ±0.1, in good 
agreement with the values of M* measured herein. 

We note that our random errors for M* in this pa- 
per are larger than the quoted errors in Mancone et al. 
(2010). This is a simple result of number statistics. 
While we only have seven clusters in this study, we had 
25 (22) clusters in our z — 1.24 (1.46) bins in Mancone 
et al. (2010). This leads to a larger random uncertainty 
in M* for this current work, although we have lower sys- 
tematic uncertainty for M* in this paper because the 
requirement of fixing a in Mancone et al. (2010) intro- 
duced a systematic uncertainty of <~ 0.2 mags into M* . 

We also compare to Strazzullo et al. (2010) who mea- 
sure the ii-band LF of a z = 1.39 galaxy cluster to 
M* + 4. They find a H = -l.2±°^ 5 , also consistent 
with our results. While there is a difference in passband 
between our studies, both trace rest-frame wavelengths 
redward of the 4000A break so we expect the difference 
in passband to have a minimal affect on our fitted values 
of a. 

Recent studies have found a deficit of low-luminosity 
red-sequence galaxies in high-redshift clusters (De Lucia 
et al. 2004; Rudnick et al. 2009; Lemaux et al. 2012). 
At face value this seems in contradiction with the flat 
a values found in this study as well as Strazzullo et al. 
(2010). However, neither our results nor the results from 
Strazzullo et al. (2010) are limited to red-sequence galax- 
ies, and therefore this apparent difference can simply be 
a sign that low-luminosity galaxies are in place in the 
cluster environment at these redshifts but have not yet 
finished transitioning onto the red sequence, as was sug- 
gested in Lemaux et al. (2012). 

4.2. Low- Redshift Comparison 

We compare our results to low-redshift cluster LFs to 
assess the evolution of a over a substantial fraction (^9 
Gyr) of cosmic history. To accomplish this we have com- 
piled a list of a measurements from the literature for 
a variety of clusters or cluster samples at different red- 
shifts, which we summarize in Table 4. We note that de 
Propris et al. (1998), Andreon (2001), and De Propris & 
Christlein (2009) did not present a formal error for a but 
did plot confidence regions for their fit, so we estimated 
the error on a from the plots of their confidence regions. 
Specifically, we derived the error from the full range of 
values covered by their la confidence contours. We split 
Table 4 up into two sections: studies which trace the rest- 
frame optical and studies which trace the rest-frame NIR 
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TABLE 4 
a Values From the Literature 



Reference Cluster a # Clusters Band (z) a 



Rest-Frame NIR 



de Propris et al. (1998) 


Coma 


1 


H 


0.023 


-0.78±0.3 b 


Andreon (2001) 


AC 118 


1 


Ks 


0.3 


-1.18±0.15 b 


Lin et al. (2004) 




93 


Ks 


0.043 


-0.84 ±0.02 


Jenkins et al. (2007) 


Coma 


1 


3.6/im 


0.023 


-1.25 ±0.05 


Muzzin et al. (2007) 




15 


K 


0.296 


-0.84 ±0.08 


Skelton et al. (2009) 


Norma 


1 


Ks 


0.016 


-1.26 ±0.1 


De Propris & Christlein (2009) 




10 


K 


0.07 


-0.98±0.2 b 


Mancone et al. (2010) 




35 


3.6/im 


0.37 


-0.60 ±0.2 


This Work 




7 


3.6/im 


1.35 


-0.97 ±0.14 


This Work 




7 


4.5/^m 


1.35 


-0.91 ±0.28 



Rest- Frame Optical 



Mobasher et al. (2003) Coma 1 R 0.023 

De Propris et al. (2003) • • • 60 Bj <0.11 

Chiboucas & Mateo (2006) Centaurus 1 V 0.0114 

Strazzullo et al. (2010) XMMU J2235-2557 1 H 1.39 



i -,0+0.04 
-J-J-8-0 .02 

-1.28 ±0.03 

_i 4+0.1 
LA -0.18 



-1.2 



+0.2 
0.15 



a A cluster name is given only when a single cluster is studied in the given paper. 
b Formal errors were not given but were estimated from plots of confidence regions. 



(such as this work). Star formation can be an important 
contributor to the rest-frame optical LF, and as such a 
is not necessarily directly comparable between the two 
sets of studies. 

The results in Table 4 are presented graphically in Fig- 
ure 4. In this Figure the fitted a values and errors are 
plotted for all the rest-frame NIR results in Table 4. 
There is substantial study-to-study scatter at low red- 
shift, and large error bars at high redshift, but Figure 4 
shows no obvious evidence for evolution in a from z = 
to z ~ 1.4, representing nearly 70% of cosmic history. 

Past work on cluster LFs have primarily characterized 
the evolution of M*, and shallow imaging has required 
assuming a value for a and fixing it as a function of red- 
shift (see, e.g., Muzzin et al. 2008 and Mancone et al. 
2010). Fixing a has been a potential source of system- 
atic uncertainty, as the strong coupling between M* and 
a means that if a is improperly held fixed then the fitted 
values of M* will also be wrong. This can be particularly 
important for studies of the evolution of M* because if 
a is evolving but assumed to be fixed then this false 
assumption can create spurious evolution in M*. This 
potential source of systematic uncertainty was discussed 
in detail in Mancone et al. (2010) because we found that 
for z > 1.3 the fitted values of M* to the cluster LF 
deviated strongly from passive evolution. We concluded 
that while evolution in a could contribute to the mea- 
sured deviations from passive evolution, it was unlikely 
to be the underlying cause because the direction of the 
deviation would require a to become steeper at higher 
redshift. Having now measured a at high redshift we 
can conclude that evolution in a was not the cause of 
our observed deviations from passive evolution as a does 
not evolve significantly with redshift out to z ~ 1.4. 

4.3. Comparison to the Field LF 

We find that at high redshift the faint-end slope of the 
cluster LF matches that of the field. Saracco et al. (2006) 



Age of Universe (Gyrs) 
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Fig. 4. — Best-fitting values of a to the cluster luminosity func- 
tion versus redshift. Data points come from the rest-frame NIR 
studies from Table 4 and include a variety of literature results. 
The results from this work are shown as red circles while previous 
results are shown as black triangles. The dashed horizontal line 
denotes a = —1.0, corresponding to a flat LF for faint galaxies. 

measured a in the rest-frame J-band for field galaxies, 
finding a = — 0.941q '\f in a redshift bin centered at z ~ 
1.2. Cirasuolo et al. (2007) found a = -0.92 ± 0.18 for 
galaxies with 1.25 < z < 1.5 in the rest- frame if-band, 
again in good agreement with the faint-end slope of the 
cluster LF found herein. Moreover, work in the field 
(Saracco et al. 2006; Cirasuolo et al. 2007; Stefanon & 
Marchesini 2011) has also found that the faint end of the 
LF for field galaxies is consistent with being constant out 
to the highest redshifts studied. For example, Kochanek 
et al. (2001) find a = -1.09 ± 0.06 for 2MASS galaxies 
at low redshift, while Stefanon & Marchesini (2011) find 
that the faint-end slope of their rest-frame J-band LF is 
consistent with flat from 1.5 < z < 3.5. They compare 
their results to lower redshift studies, finding no evidence 
for evolution in a out to z = 3.5 with a mean value of 
a = —1.05 ± 0.03. This is all consistent with our finding 
that a does not substantially evolve but is consistent 
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with flat from z = to z~l. 4 for cluster galaxies. 

4.4. Implications for Galaxy Formation 

As discussed above, we find no statistically significant 
evidence for evolution in the faint-end slope of the NIR 
luminosity function out to z ~ 1.4. A lack of evolution 
in a combined with a lack of evolution in M* (Mancone 
et al. 2010) means that both faint and bright galaxies 
are largely in place at high redshift. This places a strong 
constraint on the luminosity evolution of cluster galaxies 
at z < 1.4. Either little evolution is happening at lower 
rcdshifts, or the processes responsible for LF evolution 
have no net impact on the cluster population. 

To understand the implications of this for galaxy evo- 
lution, we must understand how the various process that 
cause cluster galaxy evolution would affect the luminos- 
ity evolution of the cluster galaxy population. The fact 
that the NIR traces old stellar populations means that 
our results are most sensitive to processes which would 
cause evolution in the stellar mass of cluster galaxies. 
Any process which affects low and high stellar mass clus- 
ter galaxies equally will lead to evolution in M* . Con- 
versely, any process which leads to differential evolution 
between galaxies with low and high stellar masses will 
lead to evolution in a. 

The two primary processes by which galaxies can grow 
their stellar masses over time are star formation and 
mergers. The downsizing paradigm (see Fontanot et al. 
2009 and references therein) suggests that star formation 
will be preferentially found in lower mass galaxies at low 
redshift. Such a mass dependence for galaxy star for- 
mation histories will necessarily imply evolution in the 
faint-end slope of the LF. The amplitude of this effect 
however depends upon the total amount of ongoing star 
formation, which will vary between clusters and may de- 
pend upon the total mass of the host cluster halo. There 
is evidence that substantial star formation is still ongo- 
ing in our cluster sample (Snyder et al. 2012), as well 
as other clusters at similar redshifts (Hilton et al. 2010; 
Tran et al. 2010; Fassbender et al. 2011). In contrast, 
Muzzin et al. (2012) find that star formation has already 
been strongly quenched in their cluster at z — 1.2. 

Mergers can also build up the stellar mass of cluster 
galaxies, although mergers are expected to be suppressed 
in the cluster environment due to the high relative veloc- 
ities of cluster galaxies (Alonso et al. 2012). Recent the- 
oretical studies (Murante et al. 2007; Conroy et al. 2007; 
Puchwein et al. 2010) suggest that for massive galax- 
ies growth by mergers becomes very inefficient (but see 
Rudnick et al. 2012), and it is possible that mergers can 
yield a steepening of the faint-end slope if this efficiency 
is strongly mass-dependent. 

In contrast gravitational interactions such as galaxy- 
galaxy interactions, interactions of a galaxy with the 
cluster potential, galaxy harassment, and the dissolution 
of cluster galaxies, can all strip mass away from low- 
mass galaxies in particular (Moore et al. 1996; Boselli & 
Gavazzi 2006; Murante et al. 2007) and therefore cause 
a to grow flatter or turn over with time. 

Another process for consideration is the infall of new 
galaxies into the cluster, the effect of which depends on 
the shape of the LF for the infalling galaxy population. 
Since we find that the cluster and field galaxy popula- 
tions have a similar faint end slope (Section 4.3), we ex- 



pect that the infall of new galaxies into the cluster will 
primarily act to mitigate any potential evolution of the 
cluster LF by driving the cluster LF back towards a flat 
faint-end slope. Instead, the infall of new galaxies will 
lead to an increase in $*, the normalization of the LF, 
which we have not constrained 

Finally, a mass-dependent galactic initial mass func- 
tion (IMF) can cause the shape of the LF to change 
relative to the underlying stellar mass function. This 
is because the rate of luminosity evolution for a stel- 
lar population depends sensitively on the IMF (Conroy 
et al. 2009). Recent work (Cappellari et al. 2012) has 
suggested that the IMF does indeed depend on galaxy 
mass, such that lower mass galaxies have a flatter IMF 
(i.e., a higher fraction of high- mass stars). A flatter IMF 
leads to a faster fading of the underlying stellar popula- 
tion (Conroy et al. 2009) and therefore, if true, the results 
of Cappellari et al. (2012) suggest that low-mass galax- 
ies should fade faster than high mass galaxies when the 
stellar masses of both remains fixed. This will cause a 
to grow flatter or turn over with time, effectively act- 
ing against processes which build up the stellar mass 
of galaxies but without impacting the underlying mass 
function. 

Clearly, there are many processes which could poten- 
tially cause evolution of the NIR cluster LF at z < 1.5. 
Therefore, the lack of evolution in a observed herein, 
combined with the lack of evolution in M* observed over 
a similar redshift range (Mancone et al. 2010), places an 
important constraint on these processes. In net, they 
cannot cause any large evolution in the shape of the NIR 
cluster LF. This could be because both low mass and 
high mass cluster galaxies are largely assembled at high 
redshift, or because the differing effects of these processes 
causes little evolution in net. The latter might imply an 
uncomfortable degree of fine tuning. In general though, 
the ability of infalling galaxies to dilute any evolution in 
the cluster LF allows for more flexibility in the strength 
of other processes. 

5. CONCLUSIONS 

We measure the 3.6 and 4.5/itm luminosity functions 
of seven galaxy clusters at 1 < z < 1.5, specifically in- 
vestigating the shape of the LF for faint galaxies. We 
find the LFs to be well-fit by Schechter distributions 
with faint-end slopes of 0:3. 6 Mm = —0.97 ± 0.14 and 
ay. 5^m = —0.91 ± 0.28, both consistent with having flat 
faint-end slopes within la. Our primary conclusions are 
summarized here: 

1. We compare to studies of the NIR LF of low- 
rcdshift clusters and find no statistically significant 
evidence for evolution of the faint-end slope of the 
cluster LF. Therefore we conclude that the faint 
end of the cluster LF has not evolved significantly 
over 70% of cosmic history. 

2. Having measured a non-evolving faint-end slope we 
have removed one source of systematic uncertainty 
from studies of the evolution of M* as a function 
of redshift. This is particularly relevant for our re- 
cent detection of deviations from passive evolution 
at high redshift (Mancone et al. 2010). Shallow 
imaging in Mancone et al. (2010) necessitated fix- 
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ing a, which could have lead to spurious evolution 
in M* if a was evolving. 

3. We compare to the faint end slope for field galaxies 
at similar redshifts and find good agreement. Field 
studies (Saracco et al. 2006; Cirasuolo et al. 2007; 
Stcfanon & Marchcsini 2011) find a faint-end slope 
consistent with flat at high redshift, and the most 
recent results (Stefanon & Marchesini 2011) find 
no evidence for evolution out to z = 3.5. 



4. Given recent studies (Muzzin et al. 2008; Mancone 
et al. 2010) which have found that the evolution 
of the bright end of the cluster LF is consistent 
with passive evolution out to z ~ 1.3, we con- 
clude that the shape of the cluster LF has been in 
place and evolved little since z ~ 1.3. This could 
suggest that low-mass galaxies are largely assem- 
bled at high redshift. Conversely, it could simply 
mean that the many processes which cause evolu- 
tion of the cluster galaxy population have no net 
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